Astronomy & Astrophysics manuscript no. aa2006'7017 


© ESO 2008 


February 5, 2008 





o 
o 

(N 
X> 

IX, 



(N 
> 

On 
00 

o 
o 

Oh 
6 



Relation between photospheric magnetic field and chromospheric 

emission 

R. Rezaei 1 , R. Schlichenmaier 1 , C. A.R. Beck 1 - 2 , J.H.M.J. Bruls 1 , and W. Schmidt 1 

1 Kiepenheuer-Institut fur Sonnenphysik, Schoneckstr. 6, 79 104 Freiburg, Germany 

2 Instituto de Astroffsica de Canarias (IAC), E 38 205, La Laguna, Espain 

Received 22 December 2006/Accepted 25 January 2007 

ABSTRACT 

Aims. We investigate the relationship between the photospheric magnetic field and the emission of the mid chromosphere of the Sun. 
Methods. We simultaneously observed the Stokes parameters of the photospheric iron line pair at 630.2 nm and the intensity profile 
of the chromospheric CanH line at 396.8 nm in a quiet Sun region at a heliocentric angle of 53°. Various line parameters have 
been deduced from the CanH line profile. The photospheric magnetic field vector has been reconstructed from an inversion of the 
measured Stokes profiles. After alignment of the Ca and Fe maps, a common mask has been created to define network and inter- 
network regions. We perform a statistical analysis of network and inter-network properties. The H-index is the integrated emission 
in a 0.1 nm band around the Ca core. We separate a non-magnetically, Hnon. and a magnetically, Hmag. heated component from a 
non-heated component, Hco in the H-index. 

Results. The average network and inter-network H-indices are equal to 12 and 10 pm, respectively. The emission in the network is 
correlated with the magnetic flux density, approaching a value of H w 10 pm for vanishing flux. The inter- network magnetic field is 
dominated by weak field strengths with values down to 200 G and has a mean absolute flux density of about 1 1 Mx cm~ 2 . 
Conclusions. We find that a dominant fraction of the calcium emission caused by the heated atmosphere in the magnetic network 
has non-magnetic origin (H m ag ~ 2 pm, H n on ~ 3 pm). Considering the effect of stray light, the contribution from an atmosphere 
with no temperature rise to the H-index (Hco ~ 6 pm) is about half of the observed H-index in the inter-network. The H-index in the 
inter-network is not correlated to any property of the photospheric magnetic field, suggesting that magnetic flux concentrations have 
a negligible role in the chromospheric heating in this region. The height range of the thermal coupling between the photosphere and 
low/mid chromosphere increases in presence of magnetic field. In addition, we demonstrate that a poor signal-to-noise level in the 
Stokes profiles leads to a significant over-estimation of the magnetic field strength. 
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1. Introduction either the calcium intensity profile (e.g.. ICram & Darnel 1 1983b 

_,, , , „ ILites et aTl 1 19931) or combinations of calc ium filtergrams and 

The dominant pattern covering the entire solar surface, except „, t _ / 15 e T tiJbnnih c ™ u 

F , . , 6 . „ magnetograms (e.g., Berger & Title 2001). From the observa- 

sunspots, is granulation the top of small-scale convection cells tion of ^ caldum tmm alon ^ it is not ible to distin . 

£ . with diameters of 1 -2 Mm. On a larger scale of ~ 20 Mm super- ifjh between ^ m ic and non netic heati co 

granules are observed which show the same pattern as granu- „ Combini calcium filtergrams with magnetograms al- 

lation: a horizontal flow from the center towards the boundary , . . .. . , . ° , • e 

, . „ , lows to separate those components, but the spectral informa- 

of the cell These large scale convection cells have a mean life- don h ^ simultaneous observations that allow to reconstruct 

time of 20 h much longer than the granular time-scale (some ^ m ic fi eld and re cord the trum for the Ca n jj line 

10 minutes). The chromospheric network forms at the boundary ^ mm ( LiteTeTd] fl999l) and only available at lower 

of the supergranular cells. The network is presumably formed ^ resolution. Therefore, it is not surprising that none of 

by the long-time adve ction of magnetic flux to the boundarie s the t theories> me chanical an d Joule heating, wa s con- 

of the supergranules JPriest et alj I2Q02J; | Cattaneo et alj ^> firmed or reje c ted observationally (IFossum & C arlsson 2005; 

The interior of the network cells, the inter- network, has muc h |^^^|200i . 
less magnetic held than the network (e.g., Kell er et al.M 19941) . 

The magnetic and thermodyna mic propertie s of the network and |Schrijver| (|1987|) studied a sample of late-type stars and in- 

inter-network are different (e.g.. lLitesll2002l) . troduced the concept of the basal flux to separate the non- 

The chromospheric heat ing mechanism is one of the m ain magnetic heating from the magnetic one. The Sun is not 

challenges of solar physics dNarain & Ulmschneider|[T996l and a ver y acti ve star, with a chromospheric radiativ e loss on 

references therein). The core emission of the Ca n H and K lines the or d er of the basal flux for Sun-like stars dFawzy et alj 

is an important source of radiative losses in the chromosphere. |2002b|). |Carlsson & StgM (|1997[) presented a semi-empirical 

Moreover, this emission is an important tool to study the tem- hydrodynamic model where enhanced chromospheric emission 

perature stratification and the magnetic acti vity of the outer at- ^ due to outward propagating acoustic wav es. This model 

mosp here of the Sun and other stars (e.g., ISchrfiver & Zwaanl was criticized in subsequent inve stigations (e.g.,|Kalkofen et alj 



2000). Most of the observational studies based on these lines use Ll999| possum & Carlsson||2005|), which was partly due to dis 

agreemen ts on the temperature stratification i n h igher lay 



Send offprint requests to: rrezaei@kis.uni-freiburg.de 



ers (e.g., lAvresi [2002; Wedemev er-Bohm et al.l 120051) . In ad- 
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dition, it is not clear whether high- or low-frequency acous- 
tic waves play th e dominant role in the energy transpor t 
to higher layers dFawzv et al.l l2002at iJefferies et alj |2006 ). 
While there are theoretical indications that the magnetic fill- 
ing factor discriminates between different regimes of heating 
( Sola nki & Steinerll 990). there is no canonic model reproducing 
thermal , dynamical and magnetic propertie s of the solar chromo- 
sphere (lJudge & PeteJll99alRuTtenlll999l) . 

The PO l arimet ri c LIttrow Spec trograph (POLIS, 
ISchmidt et ail 120031; iBeck et al.l l2005al) was designed to 
provide co-temporal and co-spatial measurements of the 
magnetic field in the photosphere and the CanH intensity 
profile. We use POLIS to address the question of the chro- 
mospheric heating mechanism, by comparing properties of 
network and inter-network in photosphere and chromosphere by 
a statistical analysis. With the information on the photospheric 
fields, we separate the contributions of the magnetically and 
non-magnetically heated component. For the first time we 
study the correlation of the chromospheric emission with the 
corresponding amplitude/area asymmetry and Stokes-V zero- 
crossing velocity at the corresponding photospheric position. 
We also investigate the magnetic field strength distribution 
of the inter-network to check whether it consists of weak 



Table 1. The definition of the characteristic parameters of the 
Ca ii H profile for the peak sample (upper part) and the band sam- 
ple (lower part) (see also Fig. 2). Wavelengths are in nm. 



fields (e.g.. iFaurobert et al.ll200lt IColladosl|200ll) or kilo-Gauss 
fields (e.g. Jsanchez Almeida et al.ll2003allbh ~ 

Observations and data reduction are explained in Sects. 2 and 

3. Histograms of the obtained parameters are presented in Sect. 

4. Correlations between the chromospheric and photospheric pa- 
rameters are addressed in Sect. 5. The heating contributions are 
elaborated in Sect. 6. Discussion and conclusions are presented 
in Sects. 7 and 8, respectively. Details of the magnetic field pa- 
rameters in the inter-network are discussed in Appendix A. The 
straylight contamination and calibration uncertainties for the 
POLIS Ca channel are estimated in A ppendix B. An ov erview 
of some of our results also appears in lRezaei et al. I d2007h . 



2. Observations 

We observed a network area close to the active region NOAA 
10675 on September 27, 2004, with the POLIS instrument at 
the German Vacuum Tower Telescope (VTT) in Tenerife. The 
observations were part of a coordinated observation campaign 
in the International Time Program, where also the Swedish 
Solar Telescope (SST), Dutch Open Telescope (DOT) and the 
Telescope Heliographique pour l'Etude du Magnetisme et des 
Instabilites Solaires (THEMIS) participated. Here we analyze a 
series of thirteen maps of the same network region at a helio- 
centric angle of 53°, taken during several hours. We achieved a 
spatial resolution of around 1 .0 arcsec, estimated from the spa- 
tial power spectrum of the intensity maps. 

All maps were recorded with a slit width of 0.5 arcsec, a slit 
height of 47.5 arcsec, and an exposure time per slit position of 
4.92 s. The scan extension for the first three maps was 40.5, 55.5, 
and 20.5 arcsec, while it was 25.5 arcsec for the remaining ten 
maps. 

The CanH line and the visible neutral iron lines at 
630.15 nm, 630.25 nm, and Tii 630.38 nm were observed with 
the blue (396.8 nm) and red (630 nm) channels of POLIS. The 
spatial sampling along the slit (y-axis in Fig. 1) was 0.29 arc- 
sec. The scanning step was 0.5 arcsec for all maps. The spectral 
sampling of 1 .92 pm for the blue channel and 1 .49 pm for the red 
channel leads to a velocity dispersion of 1.45 and 0.7kms _1 per 
pixel, respectively. The spectropolarimetric data of the red chan- 



quantity: peak sample 


description 


H 3 


core intensity 


H 2 v 


violet emission peak 


H 2 r 


red emission peak 


V/R 


H 2V /H 2r 


emission strength 


H2V/H3 


A(H 3 ) 


calcium core wavelength 


^(H 2V ) 


H 2 y wavelength 


^(H 2r ) 


H 2 r wavelength 


quantity: band sample 


description 


H-index 


396.849 + 0.050 


H 3 


396.849 + 0.008 


H 2V 


396.833 ±0.008 


H 2r 


396.865 ±0.008 


Wl 


outer wing: 396.632 ± 0.005 


W2 


middle wing: 396.713 ± 0.010 


W3 


inner wing: 396.774 ± 0.010 



nel were corrected for instrumental effe cts and telescope po lar- 
ization with the procedures described bv lBeck et al. (200 5allbl) . 

Figure 1 displays an overview of one of the thirteen maps 
after spatial alignment. The map (a) shows the Fe 1 630 nm con- 
tinuum intensity normalized to the average quiet Sun intensity. 
The map (b) shows the CanH wing intensity, taken close to 
396.490 nm (for a definition of line parameters see Fig. 2 and 
Table 1). These two maps were used for the spatial alignment of 
the red and blue channels. The next two maps (c and d) show 
the intensities in the outer and inner wings (Wl, respectively, 
W3). The inner wing samples a wavelength band close to the 
core; hence, it is more influenced by the line-core emission and 
shows higher contrast of the network than the outer wing. The 
next map (e) is the H-index, i.e., the intensity of the calcium 
core integrated over 0.1 nm (cf. Table 1). The network features 
appear broadest and show the highest contrast in this map. The 
map (f) demonstrates network and inter-network masks (Sect. 
3.1). The map (g) shows the magnetic flux density obtained from 
the spectro-polarimetric data (Sect. 3.2). 



3. Data analysis 

In this section we discuss characteristic parameters of the cal- 
cium profile and the Fe 1 630 nm line pair. We briefly explain the 
inversion method that was used to infer the vector magnetic field. 



3. 1 . Definition of network and inter-network 

For each map, we created a mask to distinguish between the net- 
work and inter-network regions (Fig. If). This was done manu- 
ally on the basis of the magnetic flux and the H-index. We did 
not use the black region in Fig. If which separates network from 
inter-network. To study structures in these two regions, we de- 
fine two statistical samples: 

The peak sample: Each member of this sample has a reg- 
ular V profile and a double reversal in the calcium core. For 
this sample, we study quantities related to the Ca 11 H emission 
peaks (Table 1 , upper part). The number of data points in the net- 
work and inter-network are 19509 and 1712, respectively. We 
emphasize that the few inter-network points in this sample are 
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Fig. 1. From left to right: a) the F e 1 630 nm continuum intensity, b) the calcium wing intensity at 396.490 nm which was calibrated 
to FTS data ( Sten floet al.l 1 1 9841) . c) the outer calcium wing intensity (Wl), d) the inner calcium wing intensity (W3), e) the H- 
index, f) the masks which separate network (white) from the inter-network (gray), and g) the magnetic flux density obtained from 
the inversion. We did not use the black region between the network and inter-network. Each small tickmark is 1 arcsec. Note that 
sampling in x and y directions are different. 




396.4 396.5 396.6 396.7 396.8 396.9 
wavelength (nm) 

Fig. 2. Sample averaged calcium profile of one of the maps (the 
average profile is similar for all other maps). The bands are ex- 
plained in Table 1 . 



not distributed uniformly. They usually belong to small-scale 
magnetic elements surrounded by large areas without magnetic 
signal above the noise level. 



The same normalization coefficient was applied to all profiles of 
the map. Hence, all profiles are normalized to the average con- 
tinuum intensity. 

Table 1 lists the characteristic parameters we define for each 
calcium profile: 

a) The intensity of the CanH core, H3, and the amplitudes of 
the emission peaks, H2V and H2r, and their respective posi- 
tions, Au 3 ,An 2v , and A^ t (peak sample). 

b) Intensities in the outer (Wl), middle (W2), and inner line 
wings (W3), which are integrated over spectral bands with a 
fixed spectral range (band sample). 

From these quantities we derive: 

1. the ratio of the emission peaks, V/R = Fhv/thr, and 

2. the emission strength which is the intensity of the violet 
emission peak divided by the core, H2V/H3. 

The band intensities are defined such that the wavelength in- 
tervals are fixed for all the profiles and do not depend on the 
peak/core positions (Table 1, lower part). Figure 3 shows a com- 
parison between the H2V parameters derived from the band and 
peak definitions. The correlation is very strong (also for other 
peak parameters). Thus, it is justified to use the values from the 
band sample to improve the statistics for the inter-network. 



The band sample: All network and inter-network points, 
marked in the masks, are present in this sample. Since some 
of these profiles do not have two emission peaks in the cal- 
cium core, we use integra ted intensities in a fixed spectral 
range (ICram & Damel[T983l) at the core and wing of the calcium 
profile (Table 1, lower part). This definition retrieves reasonable 
chromospheric parameters. In total from all the maps, 24225 
points were selected as network (white in Fig. If) and 20855 
as inter-network positions (gray in Fig. If). 

3.2. Analysis of Ca 11 H profiles 

We averaged profiles over a large area (including network) to ob- 
tain an average profile for each map. These average profiles are 
the mean of more than three thousand profiles each time; an ex- 
ample is shown in Fig. 2. We then normalize d the intensity at th e 
line wing at 396.490 nm to the FTS profile dStenflo et al.l ll984). 



3.3. Analysis of Fe 1 630 nm line profiles 

The data of the red channel of POLIS were treated with 
the standard polarim etric calibration procedures described in 
Beck et all d2005albl) . The spectral line curvature was removed 
using the tellur i c O2 li ne at 630.20 nm with a routine described 
in iRezaei et al.l (120061) . The polarization signal in Q(A), U(A), 
and V(A) is normalized to the local continuum intensity, I c , for 
each pixel. The rms noise level of the Stokes parameters in 
the continuum was <x = 8.0 xl0~ 4 I c . Only pixels with V signals 
greater than 3 <x were included in the profile analysis. For regu- 
lar Stokes-V signals above the threshold, we derived positions 
and amplitudes of the profile extrema in all Stokes parameters 
by fitting a parabola to each lobe. 

An inversion was perf o rmed using the SIR code 
dRuiz Cobq & del Toro Iniestal Il992l). We used the sam e 
setup as in lBellot Rubio & Beckl d2005h and iBeck et all d2006l) : 
a two-component solar atmosphere model with one magnetic 
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middle wing (pm) outer wing (pm) Cawing Fe continuum 

Fig. 4. Histograms of the intensity parameters for the network (thick) and inter-network (thin), using the band sample: a) the H- 
index, b) H2V, c) F^r, d) W3 (the inner wing intensity), e) W2 (the middle wing intensity), f) Wl (the outer wing intensity), g) the 
calcium wing intensity at 396.490 nm which was used for normalization (see Fig. 2), and h) Fei 630.2 nm normalized continuum 
intensity (see Table 1 for definitions). 
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0.05 0.10 0.15 0.20 0.25 

H 2v peak 

Fig. 3. Scatter plot of the H2V band vs. peak definition. In the 
band definition, a fixed spectral range was used. Gray and black 
show the network and inter-network, respectively. 

and one field-free component. Additionally, a variable amount 
of straylight was allowed for. This accounts for unresolved mag- 
netic fields inside each pixel. We did not consider any gradient 
for the atmospheric parameters except for the temperature. The 
inversion yields the magnetic field vector and an estimate for 
the magnetic filling fraction. The flux density map (magnetic 
flux per pixel, Fig. lg) is based on the inversion results. 

4. Histograms 

4.1. Parameters of the intensity profiles 

The network patches are brighter than the inter-network at all 
wavelengths: in the core and wing of the CanH line and in 
Fe 1 630 nm continuum. The distributions of the intensity param- 
eters at different wavelengths are shown in Fig. 4. Thick and thin 
curves show network and inter-network, respectively. A com- 
parison between the network and inter-network histograms in- 
dicates that the calcium core parameters (H-index, H2V, H2r) 



are more shifted than the wing/continuum intensities (Fig. 4a, 
b, c). In other words, the intensity contrast increases from the 
continuum toward the wing and peaks at the core. Moreover, 
the wing intensities show different behaviors due to their re- 
spective distances to the core: the closer the selected wave- 
length band to the core, the larger the shift between network 
and inter-network distributions. Hence, the histogram of the in- 
ner wing intensity (Fig. 4d) has a larger shift than the middle 
and outer wings (Fig. 4e and f). The distribution of the CanH 
win^Q and Fe 1 630 nm continuum intensities also show a shift 
between the network and inter-network: the peak of the distribu- 
tion in the network is slightly brighter than in the inter-network 
in 396.490 nm and also in the 630 nm continuum (Fig. 4g and h). 

4.2. Magnetic field parameters 

From the inversion, we obtain the magnetic field and the mag- 
netic flux in the photosphere. The statistics of these quantities 
in the network and inter-network are discussed here; the rela- 
tion with the chromospheric parameters is discussed in the next 
section. We restrict the analysis to those profiles with a regu- 
lar Stokes-V and a double emission calcium profile, although 
there are only 1712 profiles in the inter-network. The histogram 
of the absolute magnetic flux density, <$> = f B cos a, is shown 
in Fig. 5a, where / is the magnetic filling factor, B is the field 
strength, and a is the angle between the magnetic field orienta- 
tion and the line-of-sight. The histogram for the inter-network 
has a narrow distribution, which peaks at some ~4Mxcm~ 2 
whereas there is a long tail up to ~ 250 Mx ctrT 2 for the network. 
The distribution function of the magnetic field strength (Fig. 5b) 
has a peak around ~ 1.3 kG for the network elements and in- 
creases with decreasing field strength for the inter-network. The 
peak at 200 G is due to the detection limit of the polariza- 
tion signal. The histogram of the V amplitude (Fig. 5c) shows 
weak inter-network signals with a V amplitude far less than 
1%. The distribution of the inter-network V amplitude^ peaks 



1 Here, wing means the wavelength that was used for the intensity 
calibration (see Fig. 2). 

2 Histogram of the V-amplitude of the inter-network (Fig. 5c) in- 
cludes all data points with a clear V signal. 
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Fig. 5. Histograms of the magnetic field parameters for the network (thick) and inter-network (thin), using the peak sample: a) the 
absolute magnetic flux density, b) the field strength, c) Stokes V amplitude, d) Fe 1 630.25 nm V velocity, e) the amplitude asymmetry, 
and f) the area asymmetry. 



Table 2. Correlation coefficients between the band intensities 
and the H-index in the network and inter-network using the band 
sample. 



parameter 


network 


inter-network 


H 3 


0.88 


0.90 


H2V 


0.92 


0.94 


H2r 


0.88 


0.93 


W3 


0.77 


0.68 


W2 


0.73 


0.53 


Wl 


0.58 


0.36 



at ~ 2 X 10 I c whereas the network shows a broad distribution 
up to ~ 0.03 I c . 

The V velocity distribution is shown in Fig. 5d. Although it 
has a peak around zero both for the network and inter-network, 
the inter-network shows a larger fraction of high-velocity V pro- 
files. The amplitude and area asymmetries in the network peak 
at small positive values of 10% and 3 %, respectively (Fig. 5e 
and f). In contrast, there is no tendency to positive or negative 
values for the inter-network asymmetries. This i s in agreement 
with other studies of quiet Sun magnetic fi elds ( Sig warth et al.l 
1 1 999t ISigwarthll200 It iKhomenko et al.ll2003l) . 

5. Correlations 

5.1. Correlations between chromospheric quantities 

There is a strong correlation between the H-index, H3, H2V, H2r, 
and the inner wing intensity (with a correlation coefficient > 0.8, 
see Fig. 6). Therefore, we adopt the H-index as a proxy for the 
calcium core emission. It also has the advantage to be the cal- 
cium intensity parameter closest to the widely used filtergramsB 



T he (V/R) ratio is an indica t or for the H3 line-c ore posi- 
tion (iRutten & Uitenbroekl[T99l ICram & Darnel 1 1983I) . There 
is a strong correlation between the (V/R) ratio and the 
redshifted and blueshifted calcium profiles correspond to (V/R) 
ratios smaller and larger than one, respectively (Fig. 6, lower 
right panelfl- We find a similar correlation between the emis- 
sion strength and the (V/R) ratio: the larger the (V/R) ratio, the 
larger the emission strength. The lower left panel of Fig. 6 shows 
the correlation between the outer wing intensity and the H-index. 
While there is a correlation in the network (gray), there is no sig- 
nificant correlation in the inter-network (Table 2). We calculated 
the correlation coefficient between the H-index and all intensi- 
ties in the line wing to investigate this difference in more detail 
(Table 2). In the network, there is a stronger correlation than in 
the inter-network for all the wing bands considered. We return 
to this point in Sect. 7. 



5.2. Correlations between photospheric quantities 

There is a strong correlation between the amplitude and area 
asymmetries, both in the network and inter-network (Fig. 7, top 
left panel). In contrast, there is no correlation between either 
amplitude or area asymmetries and the V velocity (bottom left 
panel, Fig. 7). T hese are typical pro perties of the quiet Sun mag- 
netic field (e.g. JSigwarth et aTlll999l) . The right panels of Fig. 7 
show scatter plots of the V asymmetries versus V amplitude. 
There is a tendency for high amplitude V signals to have small 
asymmetries. On the other hand, histograms of the asymmetries 
peak at a positive value. Implications of these two findings will 
be discussed in Sect. 5.3.2. 



Note that filtergrams, e.g., from DOT ide Wiin et aD l2005) also 
cover parts of the line wing while we integrate over a band with a width 
of .1+0.001 nm. 



4 The calcium core position is expressed in velocity units for com- 
parison. This does not mean that it corresponds to a Doppler shift. 
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10 12 14 16 1 
H-index (pm) 



4-2 2 4 
Hj core shift (km/s) 



Fig. 6. Upper panels: correlation between the H^v and H21- based 
on the peak sample and the H-index. Lower left: correlation be- 
tween Wi and the H-index. Lower right: correlation of the (V/R) 
ratio with the calcium core position. Gray and black show the 
network and inter-network, respectively (for abbreviations, see 
Table 1). 
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Fig. 7. Top left: scatter plot of the amplitude vs. area asymme- 
tries. Bottom left: scatter plot of the area asymmetry vs. the V 
velocity. Right panels: scatter plots of the amplitude and area 
asymmetries vs. the Fe 1 630.25 nm V amplitude. Gray and black 
show the network and inter-network, respectively. 

5.3. Correlations between photospheric and chromospheric 
quantities 

5.3.1. Calcium core emission vs. magnetic flux 

The upper panel of Fig. 8 shows the relation between the pho- 
tospheric magnetic flux and the chromospheric emission for the 
network. For lower flux values (< 100 MxctrT 2 ), there is a clear 
increase of emission with flux. However, for higher magnetic 
flux densities, the H-index increases slowly. To reproduce the 
observed relation, we utilize a power law fit to the data, 



H = a<& + c, 



(1) 



where O is the absolute magnetic flux density, H is the H-index, 
a is a constant coefficient, b is the power index, and c is the 
non-magnetic contribution (see section 5.3.2). We use a variable 
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Fig. 8. Upper panel: Correlation between the H-index and the 
absolute magnetic flux density. Gray is the original data and 
black is the binned data: each point is average of 25 points. The 
middle curve shows a fit of a power law to the original data (Eq. 
1), the other two curves give the 1 —cr error range of the fit. Lower 
panel.Theie is no correlation between the H-index and the mag- 
netic flux density in inter-network. 



lower threshold of 0, 3, . ., 20 Mx cm" 2 (cf. Table 3) for the fit and 
neglect all data points with fluxes below this level. We find that 
the power exponent, b, depends strongly on the threshold. The 
higher the threshold, the better the fit curve rese mbles a straight 
line, as first reported by Skumani ch et alj d 19751) . If we keep all 
the points, including inter-network, we obtain a value for the 
power exponent b of « 0.2. 

The calcium core emission in the inter-network does not cor- 
relate with the magnetic flux, at least with the concentrated mag- 
netic field within the Zeeman sensitivity (Fig. 8, lower panel). 
The value of the offset, c, of the fit to the network (Table 3) 
for low magnetic flux densities is around 10 pm, which is con- 
sistent with the average H-index value of all the inter-network 
profiles (10pm for the peak sample). The value of the c param- 
eter changes for different thresholds of the flux density (Table 
3). Implications of this finding for the basal flux are discussed in 
Sect. 7. 

5.3.2. Calcium core emission vs. Stokes-V velocity and area 
asymmetry 

There are different behaviors for the network and inter-network 
calcium core emission with respect to V velocity and asymme- 
tries (Fig. 9). The core emission peaks at small positive ampli- 
tude/area symmetry in the network. This is caused by strong 
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Fig. 9. Correlation between the H-index and amplitude/area asymmetry and V velocity. Plusses and squares show network and 
inter-network, respectively. The binning method is similar to Fig. 8 



Table 3. The parameters of the fit to Eq. 1 to the network data. 
The first column is the threshold for the magnetic flux den- 
sity (MxcnT 2 ). By increasing the threshold, we avoid the inter- 
network intrusions. For a comparison to the peak inter-network 
flux density, see Appendix A. 



cut 


a(pm) 




b 


c(pm) 





0.73+0.05 


0.28 


±0.01 


10.1+0.1 


3 


0.71+0.05 


0.29 


±0.01 


10.1+0.1 


5 


0.61+0.04 


0.31 


±0.01 


10.3+0.1 


10 


0.22+0.02 


0.45 


±0.02 


1 1.0+ 0.1 


20 


0.14+0.06 


0.51 


±0.06 


11.3+ 0.2 



V profiles that show small asymmetries (Fig. 7, right panels). 
The highest H-index in the network corresponds to almost zero 

V velocity. In its diagram (Fig. 9), the slope of the left branch 
(blueshift) is smaller than the right one. Moreover for the strong 
upflow or downflow in the magnetic atmosphere in the network, 
the H-index decreases significantly. The H-index in the inter- 
network does not depend on any parameter of the Stokes- V pro- 
file (Fig. 9). 

Considering the fact that the amplitude and area asymmetries 
strongly correlate with each other (top left panel, Fig. 7), a nat- 
ural consequence is that the left (right) branch of the scatter plot 
of the H-index vs. amplitude asymmetry corresponds to the left 
(right) branch of the scatter plot of the H-index vs. area asym- 
metry (Fig. 9). However, if we compare only negative or posi- 
tive branches of the scatter plots of the H-index vs. V velocity 
and asymmetries in the network, we realize that the left branch 
in the V velocity does not correspond to the similar branch in 
the asymmetry plots. The left diagrams in Fig. 10 consider only 
the profiles with a negative V velocity, while in the right pan- 
els profiles with a negative area asymmetry are shown. A pos- 
itive V velocity may thus correspond to a positive or negative 
area asymmetry. There is no relation between the amplitude/area 
asymmetries and the V signal in the inter-network (right panels, 
Fig. 7). 



6. The magnetically and non-magnetically heated 
components 

The H-index includes the Hi (the minima outside the emission 
peaks), H2, and H3 spectral regions. So its formation height ex- 
tends from the higher photosphere to the middle chromosphere. 

The Ca n H & K lines are one of the main sources of the chro- 
mospheric radiative loss. We use the H-index as a proxy for 




-60 -40 -20 20 40 60 
area asymmetry (%) 
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Fig. 10. In the left column diagrams, only the points with a neg- 
ative V velocity are plotted. In the right panels, only points with 
a negative area asymmetry are plotted. Pluses and squares show 
network and inter-network, respectively. 

the emission in the low/mid chromosphere, and assume a linear 
relation between the H-index and the chromospheric radiative 
loss. To derive the contribution of magnetic fields to the chro- 
mospheric emission, we added the H-index of all points with 
magnetic flux above a given flux threshold. The fractional H- 
index, 77, is then defined by normalizing this quantity to the total 
H-index of all points in the field of viewQ: 

,(O>O )= ^ (*>*o). 

2j H (all maps) 

Table 4 lists the fractional H-index in different magnetic flux 
thresholds. Some 20 % of the total H-index is provided by strong 
flux concentrations (<t> > 50. MxcirT 2 ) while the remaining 
80 % of that is produced by the weak magnetic field and/or field 
free regions. 

The average H-index is 12.0 and 9.8 pm in the network and 
inter-network, respectively. However, it contains some contribu- 
tions from the photosphere (outside emission peaks) and a cool 
chromosphere (without temperature rise). Therefore, we decom- 
pose the observed calcium profiles in a heated and a non-heated 
component, 

Hj = Hco + Hnon + H s q- (2) 

where Hj is the average H-index of the inter-network, Hco is 
the non-heated component, Hnon is the non-magnetically heated 



5 We use all data points in this part and do not use the mask defined 
in Sect. 3. 
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Table 4. Fractional distribution of the total H-index, 77, for dif- 
ferent magnetic flux thresholds. 



®o (Mx cm 2 ) 


# 


T]{%) 


4. 


38 905 


63.7 


10. 


26071 


44.5 


50. 


11 190 


20.4 


100. 


4933 


9.2 


150. 


1709 


3.3 


200. 


329 


0.6 



component, and H s j r is the straylight contribution. The non- 
heated component of the profile originates from the line-wing 
emission (photo sphere) and a cool chromos phere, Hco- 

We use the [Holweger & Miilleri d 19741) model atmosphere 
as a proxy for those parts of the solar atmosphere that are not 
affected by any heating process. This model is very similar to 
a theoretical radiative equilibrium model, but it was not con- 
structed to strictly satisfy this condition. We prefer this semi- 
empirical model over a theoretical radiative equilibrium model 
because it results in slightly better reproduction of the CanH 
line wings in a standard disk-center atlas. We note, however, that 
3D (M)HD simulations of the solar atmosphere commonly pro- 
duce regions in the upper photosphere and chromosphere with 
temperatures significantly lower than in the Holweger-Miiller 
model. These structures, however, are generally so small that 
their evolution is strongly influenced by their surroundings, so 
that they cannot serve as a model for those parts of the solar at- 
mosphere that are not influenced by any heating processes. That 
means that they are not suitable for determining the absolute 
minimum value of the H-index. 

We perform a consistent NLTE radiative transfer computa- 
tion in plane-parallel geometry, including effects of partial fre- 
quency redistribution, in order to estimate the non-heated com- 
ponent of the H-index. We added an isothermal hydrostatic ex- 
tension up to a height of 2000 km to the model in order to be 
able to compute the core parts of the Ca n H & K lines. This may 
seem an arbitrary choice, but considering that the computed in- 
tensity is essentially zero at the line core for models with such 
low temperatures in the upper atmosphere, its influence on the 
H-index is very small. This model atmosphe re and the obtained 
profile are similar to the COOLC model of lAvres et all (1986) 
and the corresponding profile. We find a value of Hco = 5 pm at 
disk center. At a heliocentric angle of 53°, the appropriate value 
is H C o = 5.9pm. 

We find an upper limit of ~ 12 % for the straylight con- 
tamination of the observed profiles (cf. Appendix B). Therefore, 
we decompose the average H-index in the inter-network (Hj = 
9.8 pm) to obtain the non-magnetically heated component, 
Hnon: 

Hnon = Hj - H C o - H str (3) 



Hnon = 9.8 (± 1.0) - 5.9 - 1.1 (±0.1) = 2.8 (ij;?) . 

We find an average non-magnetic heating of Hnon ~ 2.8 + 
1.1 pm. To decompose the average H-index of the network, we 
also consider the magnetically heated component, Hmag, and 
use the same method: 

Hmag = Hn - Hco _ Hnon - Hstr- (4) 



where the magnetically heated contribution is the first term in 
Eq. 1. Assuming a similar non-magnetically heated contribu- 
tion, Hnon, for the network as in the inter-network (which is 
supported by observations, e.g., ISchriiverll 19871 1 19951) . we find 
an average magnetically heated component of about Hmag ~ 
1.9 + 1.4 pm, which is consistent with the difference between 
the mean H-index in the network and inter-network. Thus, the 
non-magnetically heated component has a larger contribution in 
the observed H-index than magnetically heated one, both in the 
network and inter-network. Our error estimates in Appendix B 
show that the uncertainties of the measurements do not affect 
this result. 

7. Discussion 

7.1. Thermal coupling 

All intensity parameters, from the calcium core to the 
Fe 1 630 nm continuum, show different distributions in the net- 
work and inter-network: the network patches are brighter not 
only at all observed wavelengths of the calcium profile, but 
also in the 630 nm continuum (Figs. 7and4). As noted by 
ICram & D ame (1983), this implies that the temperature fluctu- 
ations are spatially coherent between the photosphere and the 
lower chromosphere. This coherent intensity increase is presum- 
ably due to the fact that the optical depth scale is shifted down- 
wards in the evacuated magnetic part of the atmosphere. 

The correlation between the inner wing intensity and H- 
index is strong, both in the network and inter-network (Table 
2). However, the correlation between the outer wing intensity 
and the H-index is significant only in the network (Fig. 6, lower 
left panel). This indicates that the height range of the ther- 
mal coupling between the photosphere and low/mid chromo- 
sphere increases in presence of a magnetic field. There are also 
suggest ions that this coupling e xtends to the upper chromo- 
sphere dRauscher & Marcvll2006l) . 

7.2. H-index vs. magnetic flux 

The amount of energy deposited in the network depends on 
the magnetic flux density. To reproduce the dependence of cal- 
cium core emission on the magnetic flux with a power law, we 
find that the exponent derived strongly depends on the inclu- 
sion (or exclusion) of the weakest fluxes. We obtain a power 
ex ponent, b, of ab o ut 0.3 , whi ch is smaller than the va lue given 
by ISchriiver et aT] dl989l) and lHarvev & Whitel dl999h . We as- 
cribe the difference to the higher quality of our data set: better 
spatial resolution, a lower detection limit for magnetic signals, 
an accurate estimate of the magnetic flux due to the inversion of 
vector polarimetric data instead of a magnetogram, and unlike 
these authors, we do not consider an assumption for the back- 
ground component prior to the fit. Instead, we keep it as a free 
parameter. In this way, we find a background component close 
to the mean H-index of the inter-network for the weakest flux 
concentrations. 

The rate of increase of the chromospheric emission vs. mag- 
netic flux reduces for strong magnetic flux densities. This im- 
plies that for large flux concentrations, either the filling factor or 
the magnetic field strength saturates. There are indications that 
it is the filling factor: the available space is completely filled by 
expan ding flux tubes dHammerfll987t ISaarlll996c ISchriiver et alj 
1996). In contrast, a lower limit of chromospheric emission also 
exists. In the inter-network, no relation between emission and the 
photospheric fields is found (Figs. 8 and 9). The average value of 
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the H-index in the inter-network of around 10 pm corresponds to 
the offset, c, Eq. 1 . This reflects a constant contribution to the H- 
index which is pres ent even without photo spheric magnetic flux, 
in agreement with Schriiver ( 19871 Tl995l) who argued that the 
basal flux does not depend on the magnetic activity. The basal 
flux contains two components: the non-heated (cf. Sect. 6) and 
the non-magnetically heated contributions. The non-heated con- 
tribution depends on the temperature stratification, so that we 
speculate that the non-magnetically heated component has an in- 
verse dependence on the temperature stratification. 

Figure 8 (upper panel) shows that there is a variable lower 
limit for the H-index versus the magnetic flux density. In con- 
trast, the upper limit is less clearly defined, and would be in 
agreement with a constant maximum value independent of the 
amount of magnetic flux. There are also similar behaviors for 
the upper and lower limits in the scatter plots of the H3, H2V, 
and H2r versus the magnetic flux density. This is similar to up- 
per and lower limits of the basal flux of stars vs. the color (B-V) 
where the lower boundary ch anges but the upper one is almost 
constant (Fawzv et al. 2002b, their Fig. 2). The situation for the 
inter-network is different: both the upper and lower limits are 
independent of the magnetic flux density. 



7.3. Decomposing the Ca 11 H profile 

We interpret the solar Ca 11 H line profile as the superposition of 
a cool chromosphe ric prof i le, the heated component and stray- 
light components. lOranjd (11983) first applied this method to 
reproduce diffe rent observed Ca II K profiles from two "basic" 
profiles. Later, ISolanki et al.l dl99ll) used the same idea to de- 
compose the observed profiles by a combination of two theoret- 
ically calculated profiles. Our one-dimensional NLTE radiative 
transfer calculations reveal that the H-index of a cool chromo- 
sphere is about 5.9 pm. Since the chromospheric emission of the 
Sun as a star is slightly above the minimum among the Sun- 
like stars (Sch riiver & Zw aan 2000), the average quiet Sun pro- 
file contains a heated component. Subtraction of the cool chro- 
mospheric component from the measured H-index (along with 
straylight considerations) leads to an estimation of the pure non- 
magnetically and magnetically heated components in the mea- 
sured H-index. The non-magnetically heated component is al- 
most 50 % larger than the magnetically heated component. The 
presence of a significant non-heated contribution in the H-index 
indicates that not all of the chromospheric emission emerges 
from a hot ch r omosp here, in contrast to findings from, e.g., 
Kalkofen etaD (fl999l) . The minimum heated component in a 
Ca prof ile provides evidence how coo l the chromosphere may 
be (e.g., Wedemever-Bohm et al. 2005). Moreover, it favors the- 
ories in which the non- magne tic heating plays the dominant role 
in the chromospheric heating dFawzv et alj|2002ab . 

The fact that the non-magnetic chromospheric heating con- 
tributes significantly to the chromospheric energy balance is also 
found in the relative contributions of mainly field-free inter- 
network and network areas to the total emission in the field of 
view. Magnetic flux densities above 50 Mx cirT 2 add only 20 % 
to the total emission, which is also seen in the ratio of the mean 
H-index in the network and inter-network, 12.0/9.8 « 1.2 (Table 
4). 



7.4. Uncertainties 

The main uncertainty in the non-magnetically and magnetically 
heated components is the presence of some heating in our cool 



chromosphere profile. The Holwe ger & Mullerld 19741) model at- 
mosphere is close to an atmosphere in radiative equilibrium 
without non-magnetic or magnetic heating, so no heated com- 
ponent (cf. Sect. 6) is expected to be found on its Ca profile. 
Since there i s no general agreement how cool t he chromosphere 
may be (e.g.. Kalkofen et al.lll999tlAvresll2002l) . the best answer 
would be to use t he observed calcium profile with the lowest 
H-index. However. I Solanki et al.l dl99ll) concluded that it is not 
possible to obser ve a low ac t ivity p rofile as presented by the 
COOLC model of lAvres et al.ldl986|). which is sim i lar to the one 
employed by us on base of the lHolweger & Miilled (1 19741) model 
atmosphere. Thus, our cool profile can well serve as lower limit 
of the non-heated contribution to the H-index. Note that a larger 
value for Hco would increase the contribution of the magneti- 
cally heated component, leading to a larger fractional contribu- 
tion for the magnetically heated to the total heated component. 

There may be some mixed polarity fields below our polari- 
metric detection limit which may influence the ratio of the mag- 
netically heated to the non-magnetically heated component. 



7.5. H-index vs. V asymmetries 

From the relations of the chromospheric emission to quantities 
of the photospheric magnetic field in the network, we find that 
the chromospheric emission peaks at small positive values of the 
Stokes V asymmetries and at zero Stokes V velocity (Fig. 5). It 
reflects the histograms of the respective photospheric field quan- 
tities which show sim ilar distributions even at disk center (e.g., 
ISigwarthet ail l999). However, in combination with the depen- 
dence of the emission on the flux density it is inferred that the 
stronger flux concentrations mainly show small material flows 
along with non-zero asymmetries. 

The magnetically heated component is related to Stokes 
V profiles with non-zero area and amplitude asymmetries 
(Figs. 5 and 9). A possible explanation for this finding would be 
the absorption of upward propagating acoustic waves, generated 
by the turbulent convection, by the inclined fields of expand- 
ing flux tubes. This would imply that the energy is deposited 
at the outer boundary or in the canopy of flux concentrations 
rather than in the central, more vertical, part. It is mainly because 
the non-magnetic cutoff frequency is lowered at the boundary of 
the flux tube s, where the field lines are inclined. This was first 
predicted by Suem atsul ( 19901) and recently achieved some ob - 
servational support ( Hansteen et alj|2006t iJefferies et alj |2006). 
Figure 9 indicates that the maximum observed H-index has non- 
zero V asymmetry. These asymmetric V profiles are consistent 
with the case when the line of sight passes through the canopy 
of a magnetic element or through a flux tube ax is (positive and 
negative asymm etries respectively , ISteineri [T999). Therefore our 
finding supports ISuematsu] dl990h . 



7.6. Inter-network field strength distribution 

As a byproduct of our study on the relation between photo- 
spheric fields and chromospheric emission, we obtained distri- 
butions of field strength for network and inter-network regions. 
For the inter-network fields, for the first time an inversion of vis- 
ible spectral lines in the weak field limit led to the sa me distri- 
butio n as results from the more sensitive infrared lines dColladosI 
120011) . This is discussed in Appendix A. Thus, we find that the 
inter-network magnetic field is dominated by field strengths of 
the weak field regime which ranges up to some 600 G. This 
rules out the argument that the inter-network is dominated by 
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a magnetic field with a streng th of more than a kilo-Gauss (e.g., 
ISanchez Almeida et alj2003al) . 



8. Conclusions 

There is no correlation between the H-index and magnetic field 
parameters in the inter-network. Therefore, we conclude that the 
magnetic field has a negligible role in the chromospheric heat- 
ing in the inter-network. On the other hand, the H-index is a 
power law function of the magnetic flux density in the network 
with a power index of 0.3. The average H-index observed in the 
network and inter-network are ~ 10 and 12 pm, respectively. We 
find a non-magnetic component in the network (based on Eq. 
1) of about 10 pm which shows the consistency of our analysis: 
the non-magnetic part of the network H-index is equal to the 
H-index of the inter-network. 

A NLTE radiative transfer calculation, using the Holweger- 
Miiller model atmosphere, indicates that the non-heated com- 
ponent of the H-index, emerging from a cool chromosphere, is 
about 5.9 pm. Comparison of this non-heated component and the 
average H-index in the inter-network has two implications: a) 
some of the observed chromospheric emission does not originate 
from a hot chromosphere, and b) the non-magnetically heated 
component is about 50 % larger than the magnetically heated 
component. From this, we conclude that the non-magnetically 
heated component has a larger contribution in the chromospheric 
radiative loss than the magnetically heated component, both in 
the network and inter-network. 

In our statistical ensemble, spatial positions with strong mag- 
netic field (O > 50 Mx cm 2 ) contribute about 20 % of the total 
H-index. Correlations and histograms of the different intensity 
bands in the CanH spectrum indicate that above a magnetic 
threshold, photosphere and low/mid chromosphere are thermally 
coupled. Moreover, our findings are consistent with the idea that 
the energy transfer in a flux tube has a skin effect: the energy 
transfer is more efficient in the flux tube boundary (canopy) than 
at its center (axis). 

For the first time, we find a magnetic field distribution in the 
inter-network using visible lines which is similar to results in- 
ferred from infrared lines with larger Zeeman splitting. The dis- 
tribution function increases with decreasing field strength. The 
peak at 200 G is due to the detection limit of the polarization 
signal. The average and distribution peak of the absolute flux 
density (of magnetic profiles) in the inter-network are ~ 1 1 and 
4 MxctrT 2 , respectively. We conclude that the combination of 
high spatial resolution and polarimetric accuracy is sufficient to 
reconcile the different results on field strength found from in- 
frared and visible lines. 
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Fig. A.l. Close up view of the flux density distribution for the 
inter-network. At a value of ~ 4.0 Mx crrT 2 , the histogram shows 
a peak. 
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Appendix A: The magnetic field distribution 

Traditional magnetic field strength distributions for the quiet 
Sun (inter-network) based o n visible fines have a broad 
peak centered around 1 kG ([Sanchez Almeida et al.1 l2003at 



Domihguez Cerden a et al. 2006). On the other hand, observa- 
tions in the infrared lines show a clear peak at lower field 
strengths and an almost exponentially decreasing probability for 
highe r fields strength s (Collados 2 00 it iBellot Rubio & Colladosl 
12003b iKhomenko et al.l 120031) . This discrepancy between the 
retrieved parameters fr om the visible and infrared lines has 
not b een solved so far (ISteinerl l2003t iMartfnez Gonzalez et al.l 
120061) . 

We achieved a magnetic field distribution for the inter- 
network using the Fe 1 630 nm pair that is similar to the infrared 
studies (Fig. 5.b). We elaborate on the effects of spatial res- 
olution and signal-to-noise ratio on this finding. We used the 
Kiepenheuer Adaptive Optics Syst em (KAOS) to improve s pa- 
tial resolution and image stability dvon der Liihe et aL 1 120031) . It 
provides stable sharp images in a certain field of view, and there- 
fore allows us to increase the exposure time to have better signal- 
to-noise ratio. 

Figure A.l is a close-up view of the inter-network magnetic 
flux distribution. The lowest flux detected is ~ 0.05 MxctrT 2 . 
The distribution of the flux decreases sharply for fluxes below 
cD < 4 Mx cm 2 , which is due to the detection limit of the po- 
larization signal. The small fluxes still detected give us, how- 
ever, confidence that the magnetic field strength distribution of 
Fig. 5.b is reliable as a consequence of the combination of high 
spatial resolution and good signal-to-noise ratio. 

To investigate the influence of noise on the inversion results, 
we compare inversion results of the original datasets with and 
without adding noise. In order to create "noisy" data sets, the 
rms value of the noise in the original data was increased to 
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Fig.A.2. The upper four plots show an example of the origi- 
nal Stokes profiles and the inversion fits. The lower four figures 
show the same profile with a rms noise value two times larger 
and corresponding fits. The magnetic field strength, magnetic 
filling factor, and macro turbulence velocity for the original and 
noisy profiles are (842, 1543) G, (7.6, 5.3)%, and (2.95, 2.77) 
km s _1 respectively. The x-axis shows the spectral coordinate in 
mA. 



<x= 15.0 xl0~ 4 4, two times larger than the original value. The 
spatial sampling was identical to the original data. We emphasize 
that inversion was done using the same assumptions and initial 
model atmosphere. Figure A. 2 compares two sets of full Stokes 
profiles (black) and the inversion fit (red). The upper four panels 
are the original profiles and the lower four panels are the same 
profiles with noise. The fits and retrieved model atmospheres dif- 
fer significantly. Especially important is the fact that the noisy 
inversion leads to a magnetic field of 1 .5 kG. 

To compare the field strength obtained from the original and 
noisy datasets, their histograms are shown in Fig. A. 3. The distri- 
bution of the original field strength peaks around ~ 200 G (black 
histogram), while for the noisy data, it has a clear shift toward 
higher values with a peak at 0.8 kG (red histogram). This em- 
phasizes the role of the noise in the exis ting discrepancy between 
visibl e and infrared measurements (Bell ot Rubio & Colladosl 
120031) . 

In brief, we find that higher spatial resolution of these ob- 
servations along with low noise data is an important step to- 
ward resolving disagreements between the visible and infrared 
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Fig. A.3. Histograms of the inter-network field strengths of one 
of the maps based on the inversion of original data (black) and 
noisy data (red). 



polarimetric measurements. The visible Fei630nm pair shows 
higher Stokes amplitudes for small magnetic fields than the in- 
frared lines of Fe 1 1 .56 //m. Hence, with an equal amount of 
noise in observational data, the visible lines provide a higher 
signal-to-noise ratio than infrare d lines. Therefore, in contrast to 
Mart inez Gonzalez et alj d2006l) . we find the visible Fe 1 630 nm 
pair to be a proper tool to investigate the inter-network magnetic 
field. We ascribe the present disagreements to the low spatial 
resolution of previous observations and different signal-to-noise 
ratios. 



Appendix B: Calibration uncertainties in POLIS Ca 
channel 

The POlarimetric LIttrow Spectrograph (POLIS) was designed 
to facilitate co-temporal and co-spatial measurements of the 
vector magnetic field in the photosphere and th e Ca n H intensity 
profile dSchmidt et alj|2003t iBeck et alj|2005al) . We try to elab- 
orate on the uncertainties in the intensity calibration, because it 
is critical for the estimates of the various contributions (stray- 
light, non-magnetic heating, etc.) to the observed calcium pro- 
files. Observations are discussed in Sect. B.l. We study the lin- 
earity of the CCD camera for different light levels in Sect. B.2. 
Then, we investigate instrumental and solar straylight contribu- 
tions in Sect. B.3. In Sect. B.4, we perform standard error prop- 
agation to quantify uncertainties in our final calibrated spectra. 

B.1. Observations 

In order to study linearity of the POLIS CCD camera of the cal- 
cium channel at different light levels, we observed a set of cal- 
cium spectra and fiat field data with different accumulations and 
slit widths at disk center. A sunspot was also observed with a 
similar setup and an exposure time of 0.82 s and six accumula- 
tions. Besides this, we recorded dark current data using either the 
AO or POLIS field stop. We use these two data sets to estimate 
the solar and instrumental straylight, respectively. 

All spectra and profiles in this Appendix are dark subtracted 
raw data. No other calibration process like fiat fielding was ap- 
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Fig. B.l. Comparison between two different slit widths and the 
same modulation speed. The black and blue profiles correspond 
to slit widths of 0.48 and 0.18 arcsec, respectively. The red pro- 
file is the blue profile multiplied by a constant factor of 2.6 for 
all pixels. The green profiles show the 1-cr RMS value (~ 1%) 
from the mean curve (red). 
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Fig. B.2. Dark subtracted spectrum of the umbra in a sunspot. 



plied. In order to show intrinsic noise in the data, we did not 
spatially average darks or spectra. 

B.2. Linearity 

In order to check the linearity of the calcium camera from pro- 
files by, e.g., comparing the ratio of the core to wing for differ- 
ent exposures, we recorded s pectra with different slit widths and 
identical modulation speed (Mul likin et al.l fl994l) . Figure B.l 
shows a comparison between the two flat field profiles taken on 
disk center with identical illumination, but different slit widths. 
The black and blue profiles correspond to a slit width of 0.48 
and 0.18 arcsec, respectively. If one multiplies the blue curve 
with the ratio of the slit widths of 2.67, the resulting curve agrees 
well with the profile of the broad slit for all wavelengths includ- 
ing the line core. We calculated the ratio of the two profiles with 
different slit widths as function of wavelength. The rms varia- 
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Fig. B.3. Left: a single Ca n H spectrum of the umbra and quiet sun profiles in a sunspot. Asterisks show spectral bands that were 
used to estimate the true signal (see Sect. B.3.2). Right: close up view of th e quiet sun and umb ral profiles near the calcium line 
core. The true signal (about 5 % in the wing close to the Ca core according to Wal lace et al.l d2000l) was subtracted from the umbral 
profile. The red curve shows 12 % of the quiet sun profile. The two cyan curves show an uncertainty interval of one sigma (1 %). 
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Fig. B.4. Ratio of the umbral profile to a quiet sun profile. If 
one considers all the spectrum, the ratio will be 10 + 3 % while 
considering only parts of the wing close to the Ca n H core, the 
ratio will be 12 + 1%. Note to the mismatch between the quiet 
sun and umbral profiles at spectral line positions due to different 
broadenings. 



tion of the ratio was 1%. This can be used as an upper limit for 
non-linearity effects of the CCD camera. 

B.3. Straylight 
B.3.1. False light 

We compare two sets of dark current images of the POLIS Fe 
channel: one with a field stop right after the focal plane, and 
another one with the field stop in front of the camera. In the sec- 
ond case, all light sources inside the observing room are fully 
blocked. The difference between the two measurements is about 
0.9 + 0.5 %. This gives an upper limit for the instrumental stray- 
light in the POLIS Ca channel as well, since most of the sources 
in the observing room will not emit in blue as efficient as red. 
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Fig. B.5. Distribution of the readout error in a set of 50 dark 
currents with 8 accumulations. 



B.3.2. Atmospheric and telescopic straylight 

The calcium profile in the sunspot umbra is significantly weaker 
than a quiet sun profile. However, t here is a well known s ingle 
emission peak in the calcium core (lLinskv&Avrettlll970h . We 
thus compare an inner wing band close to the calcium core, but 
excluding the core itself (asterisks in the left panel of Fig. B.3), 
in an umbral profile and a normal profile from outside the spot 
to derive an upper limit for the solar straylight contribution to 
observed profiles. 

Figure B.2 shows a slit spectrum of the umbra. To keep the 
sample profiles as uniform as possible, we chose two identical 
positions along the slit (at about 30 arcsec, Fig. B.2). The spec- 
trum of the quiet sun and the minimum spectrum of the umbra 
are shown in Fig. B.3 (le ft panel). 

Walla ce et al.l (2000) presented a Fourier transform spectral 
atlas of the umbra. In this atlas, the wing close to the um- 
bral core has an intensity of 5 % in units of the umbral contin- 
uum, while the wing close to the Fe i, ]ionTii, and Ni i lines at 
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about 396.5 nm has an intensity of 38 %. In our umbral profile, 
this continuum has an intensity of about 80 counts (asterisks, 
Fig. B.3, left panel), whereas close to the core the profile has 
about 45 counts. Assuming no straylight in the atlas profile, only 
11 counts (~22% of the observed valued) should be measured 
in the core. Thus, 78% of the observed signal have to be due to 
straylight. To reproduce 78% of the umbral profile, a straylight 
contribution of around 12 % is needed. If we assume that there is 
no real calcium signal in the umbra (which is wrong) and take the 
observed calcium profile of the umbra as pure scattered light, the 
ratio changes from 12 % to 15 %. Therefore, uncertainties about 
the scattered light in the atlas profile have minor importance. 

A close-up view of the spectral region close to the calcium 
core is shown in Fig. B.3 (right panel). For the reason mentioned 
above, in the right panel of Fig. B.3, the umbral profile was la- 
beled umbra (78%). The red curve shows 12 % of the quiet sun 
profile. The two cyan curves also show an uncertainty interval of 
one sigma (1 %). Figure B.4 shows the ratio of the umbra to the 
quiet sun profile. If we calculate mean and RMS of the whole 
spectrum, it gives a value of 10 + 3 % while considering only on 
the wing band close to the calcium core, we obtain 12+ 1 %. 
We take this value as an upper limit for the total straylight in 
the POLIS calcium channel. In Sect. B.3.1, we concluded that 
the instrumental straylight is about 1 %. Hence, most of the ob- 
tained straylight originates from the scattered solar light in the 
telescope and the earth atmosphere. 

B.4. Noise 



There are three noise sources in a CCD camera (Mullikin et al. 

119941) : 

1 . readout noise because of the pre-amplifier 

2. dark current noise which is influenced by the chip tempera- 
ture 

3. photon noise due to the quantum nature of light which obeys 
Poisson statistics. 

The readout error of the camera was estimated by a set of 
50 consecutive dark frames with similar exposure and accumu- 
lations. Here, we compare the number of counts on each pixel 
in different images. We attribute the standard deviation of the 
count value mainly to the readout error, and hence, obtain a 
readout error for each pixel separately. The distribution of the 
readout errors is shown in Fig. B.5. The readout noise is about 

^readout ~ 2 counts - 

To investigate the thermal noise in the dark current, one has 

to remove small systematic offsets between CCD columns (prob- 
ably due to thermal fluctuations of the AD system) before con- 
sidering statistical measures, e.g., rms, of a dark profile (a CCD 
row). As it is shown in Fig. B.6, for a dark frame with 6 accumu- 
lations, the resulting rms (after removing systematics) is far less 
than the standard deviation of the whole map including a large- 
scale variation. It provides an estimate for the thermal noise in 
the camera. So the rms of the dark frame is about o~thermal ~ 3 
counts. 

The total noise in the recorded data, o" no j se , consists of three 
components: 



noise 



^photon + ^thermal + ^readout 



(B.l) 



where the photon noise originates from the Poisson fluctuations 
of the number of photons counted. Assuming that np^^ is the 
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Fig. B.6. C/pper panel: a row in a dark frame with 6 accumu- 
lation. The red curve is a polynomial fit to subtract systematic 
shift from the random noise. Lower panel: The residual signal 
after subtraction of the polynomial, which is the real noise. Note 
to the large difference between the rms of the map (excluding its 
borders) and the rms of the residuals. 

n photon- ca l c i um core ls more influenced by the noise than 
the other points in the spectrum. Figure B.7 shows a dark sub- 
tracted calcium spectrum with 6 accumulations and a narrow slit 
(0.18 arcsec). The minimum number of counts in the core is 117, 
so ^photon ~ !!• Therefore, the total noise in the core (Eq. B.l) 
will be o" no j se ~ 1 1 . So the worst signal-to-noise ratio at the 
calcium core position will be better than 10. Note that the data 
described in the paper were obtained with a wide slit (0.48 arc- 
sec) while the spectrum shown in Fig. B.7 was recorded with a 
narrow slit (0.18 arcsec), so we have better signal-to-noise ratio 
in the data. 



B.5. Error propagation 

The relative intensity of the core to wing is defined as: 
Icore 



(B.2) 



wing 



number of detected photons, the corresponding rms, cr 2 



photon 



To estimate their RMS errors, we find the signal-to-noise ra- 
tio for the core and wing intensities in some wing bands. Since it 
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Fig.B.7. A sample calcium profile. We use the calcium core 
minimum counts to estimate the photon noise. In this case, the 
signal-to-noise ratio at the calcium core is more than 10. 



is not possible to use the core intensity to obtain the noise due to 
intrinsic changes, we use the average signal-to-noise of the left 
and right wings close to the calcium core as the signal-to-noise 
of the core. So the real signal-to-noise of the core is probably 
worse than that. The S/N ratio of the core and wing vary from 
35 and 112 for two accumulations to 174 and 216 for twenty 
accumulations. We use the worse S/N ratio corresponding to an 
exposure time of 1 .64 s and obtain the following values for the 
RMS errors of the core and wing intensities: 

rr„ 1 

— ■ ~ 3 % (B.3) 



Icore 34 

' -1% (B.4) 



I • 112 
*wing LLZ - 

So, we can obtain the cr c as follows: 



-= J( 7 ^) 2 + ( 7 ^-) 2 ~3%. (B.5) 
x ^ Icore l w ing 

Error sources of the final calibrated data include uncertain- 
ties from the flat-field parameter, the straylight and the linearity 
of the detector. The straylight uncertainty was obtained about 
1 % in Sect. B.3. Assuming an upper limit for the uncertainty in 
the flat-field to be 5 % and some 1 % uncertainty in linearity of 
the core to wing ratio, an upper limit for the uncertainty in the 
calibrated spectrum will be: 



cr, 

— < 10%. 
r 



